A dense micro-cluster of Class protostars in NGC 2264 D-MM1 
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We present sensitive and high angular resolution (~ 1") 1.3 mm continuum 
observations of the dusty core D-MM1 in the Spokes cluster in NGC 2264 using 
the Submillimeter Array. A dense micro-cluster of seven Class sources was de- 
tected in a 20" x 20" region with masses between 0.4 to 1.2 M and deconvolved 
sizes of about 600 AU. We interpret the 1.3 mm emission as arising from the en- 
velopes of the Class protostellar sources. The mean separation of the 11 known 
sources (SMA Class and previously known infrared sources) within D-MM1 is 
considerably smaller than the characteristic spacing between sources in the larger 
Spokes cluster and is consistent with hierarchical thermal fragmentation of the 
dense molecular gas in this region. 
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Introduction 



NGC 2264 is a very well studied hierarchically structured young embedded cluster (Lada 



fc Lada|2003 ) associated with a giant molecular cloud in the Monoceros OBI complex, 800 pc 



(±50 pc; Sung et al. 1997 Peretto et al. 2006) away from the Sun. There is abundant evi 



dence of ongoing star forming activity in this cluster, such as molecular outflows (Margulis 



et al.|1988[ [Wolf-Chase et al.|2003[ ) and Herbig-Haro objects ( |Reipurth et al.|2004[ ). |Margulis 
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et al. (1989) identified luminous far infrared sources with the InfraRed Astronomical Satellite 



(IRAS), many of them being Class and Class I objects. Submillimeter observations of these 



sources by Williams & Garland (2002), Wolf-Chase et al. (2003), and Peretto et al. (2006) 



identified several dense clumps embedded in dusty filamentary fingers of molecular material 
that radiated out from IRAS 12. More recent observations using the Spitzer Space Tele- 
scope revealed that many of these clumps contained bright 24 /im Class 0/1 sources and the 



clustering was named the Spokes cluster due to its geometrical configuration (Teixeira et al. 



2006). Using nearest neighbor analysis, Teixeira et al. (2006) found that these protostars 



have a characteristic spacing that corresponds to the Jeans length and are very likely tracing 
the primordial substructure of the cluster. They interpret the regular spacing of the proto- 
stars as a fossil signature of thermal fragmentation. The most massive and densest of these 
submillimeter cores, D-MM1 (also known as IRAS 12 SI) has a mean density of ~10 6 cm~ 3 



(Wolf-Chase et al. 2003 Peretto et al. 2006). HCO+ observations of D-MM1 by Williams & 



Garland (2002) shows a deep red-shifted absorption indicating that it is collapsing at a speed 
that Peretto et al. (2006) measured to be ~0.1kms _1 . Wolf-Chase et al. (2003) identified 



D-MM1 as a Class object harbouring one or more protostars. Recent near-infrared imag- 
ing using the PANIC camera at the Magellan telescope, as well as Spitzer IRAC imaging 
resolved D-MM1 into 11 f^-band and 8 8 fim sources ( Young et al.| [2006) within the field we 
observed with the Submillimeter Array (SMA). These high angular resolution observations 
resolved D-MM1 into a dense micro-cluster of 7 millimeter compact sources associated with 
young Class objects. We describe our observations in Section § 2 and present our results in 
Section § 3. Finally, we discuss our findings on the characterization of these new millimeter 
sources in Section 5 4. 



2. Observations 

The observations were made with the Submillimeter Array (SMA) during 2005 De- 
cember 25 in its "compact" configuration, and 2006 February 5 and 10 in its "extended" 
configuration, with eight antennas in both configurations. The frequency was centered at 
230.538 GHz in the lower sideband, while the upper sideband was centered at 220.538 GHz. 
The primary-beam size (half-power beamwidth) of the 6 m diameter antennas at 230 GHz is 
54". The phase reference centers of the field is (a, £)(J2000) = (06 h 41 m 06.45 s , +09 o 33'47.91"). 
The zenith opacity, (t23ogh z )> measured with the National Radio Astronomy Observatory 
(NRAO) tipping radiometer located at the Caltech Submillimeter Observatory, varied be- 
tween 0.03 to and 0.04 in the extended configuration observation time, and in the compact 
configuration observation time was very stable of 0.2. The phase and amplitude calibrators 
were the quasars 0739+016 and 0530+135 with measured flux densities of 2.82 ± O.Uy and 
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0.84 ± 0.1 Jy, respectively. The uncertainty in the flux scale is estimated to be 20%, based 
on the SMA monitoring of quasars. Observations of Uranus provided the absolute scale 



for the flux density calibration. Further technical descriptions of the SMA are found in Ho 



et al. (2004). The data were calibrated using the Interactive Data Language superset MIR, 



based on the MMA software developed for the Owens Valley Radio Observatory (Scoville 



et al.||l993 ) which was adapted for the SMA[] The calibrated data were imaged and analyzed 



in the standard manner using the Multichannel Image Reconstruction Image Analysis and 
Display (MIRIAD) package. We weighted the (u,v) data using the ROBUST parameter of 
INVERT MIRIAD task set to 2, optimizing for a maximum sensitivity in the continuum 
image. This option is recommended to achieve the largest signal-to-noise ratio possible, 
although some angular resolution is sacrificed. The synthesized beam had dimensions of 
l'/4 x 1'.'3 with a P.A. = 10.4°. Finally, the resulting continuum image rms noise level (a) 
was 3.5 mJy beam -1 . 



3. Results: new millimeter sources 

We detected seven compact sources whose peak fluxes exceeded 5 a within the core D- 
MM1. Figure [I] shows a panel of near-infrared (NIR) and mid-infrared (MIR) images of the 
core with the SMA sources overplotted (with contours). The i^-band image (0.125"pixel -1 ), 
depicts emission in the vicinity of the SMA sources, however this NIR emission is generally 
not pointlike (except for the source associated with SMA-3) and may correspond to scattered 



light from a central protostar. The Spitzer Infrared Array Camera (IRAC) image (Teixeira 



et al.||2006| ) (1.2"pixel _1 ), shows that SMA-1 and SMA-3 have MIR counterparts. The 
sources SMA-2, SMA-6, and SMA-7 are also associated with diffuse MIR emission. The 
point-like sources detected in the K s -band and the 8 /im images, that were not detected 



by our SMA observations, had been previously identified as Class I sources by Young et al. 



(2006) (sources # 30, 33, 35, 44, 45, 49). In Table [T] we present the physical parameters 



of these new objects. To measure the fluxes for the SMA sources, we used the NRAO 
Astronomical Image Processing System (AIPS) software. Flux densities, as well as source 
positions and sizes, were obtained using the AIPS IMFIT procedure, where each source was 
modeled with two-dimensional elliptical Gaussians. The errors cited for the flux densities 
correspond to statistical uncertainties and are dominated by the calibrational uncertainty 
that is ~ 20% as previously mentioned. 

To determine the mass we assume the emission arises from an optically thin, unresolved, 
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source that can be characterized by a single temperature (e.g. Mundy et al.||1995 Bally et al. 
1998). The mass is calculated by M = F u d 2 /(B u (T d )n u ) where v is the frequency, T d is the 



dust temperature, F v is the integrated flux density of the source, B u is the Planck function, 



and finally, k u is the dust mass opacity. We derive the latter value from Beckwith et al. (1990 
k u = 0.1(z//1200 GHz)^ cm 2 g -1 ) as being k v = (0.0192)^ cm 2 g -1 . The value of k v assumes 
a gas-to-dust ratio of 100 (Hildebrand 1983), which may not be the most adequate to use 
for protostellar sources since dust settling to the mid-plane of the disk and erosion of the 



circumstellar envelope by photo-dissociation may decrease the gas-to-dust ratio (Williams 



et al. 


2005 




Throop & Bally 


2005 


)■ 


uncertain ( 


Beckwith et al.|1990 


). The 



Additionally, the value of the emissivity index, j3, is 



from compact envelopes, hence we adopt (3 = 1.5 to calculate the mass of the sources (Wolf- 



Chase et al. 2003 Peretto et al. 2007). We use 23 K for the mass determination, as this was 



temperature measured by Wolf-Chase et al. (2003) for D-MM1. Table [T] shows the masses 
we determined, as well as the average number density, n, given by n = 3M/(47r/miH-Rg eom ), 
where /i is the mean molecular weight, 2.34, m H is the atomic H mass, and R geom is the 
geometric mean radius. Due to the uncertainties referred to above the values of the derived 
masses are good within a factor of two. The values of n should be interpreted only as an 
estimate since the error associated with the size of the source is uncertain and reflects the 
particular (u, v) coverage of the data. 



4. Discussion 
4.1. On the nature of the new micro-cluster 

Following the number counts of extragalactic objects at millimeter wavelengths from 



Maloney et al. (2005), we estimate that the expected number of 1.3mm background sources, 
above a flux density of S=20mJy, is ( N ) ~ 5.76^g^° deg -2 . The probability of finding a 
source with flux density equal or larger than 20 mJy (see Table [l|) in a 20" x 20" region is thus 
quite low, ~ 10~ 4 . We therefore conclude that all the observed SMA sources are associated 
with NGC 2264 D-MM1 core. We find a very similar value using the number counts for sub- 



millimeter galaxies reported by Laurent et al. (2005). Out of the seven millimeter continuum 



compact sources, two are associated with pointlike emission and three are associated with 
diffuse emission, (see Fig. [l]and Table [l}. Sources SMA-4 and SMA-5 do not appear to have 
any NIR or MIR emission, although there is NIR and MIR diffuse emission located between 
these sources (which could be scattered light from one or from both sources). The decon- 
volved sizes (~ 600 AU) of these objects, as well as the high densities (~10 8 ) suggest the 
1.3 mm emission is arising from compact dense regions such as circumstellar envelopes. The 
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derived masses (0.4 to 1.2 M ) are comparable to those of low mass stellar objects indicating 
that a substantial amount of the total protostellar material is in the circumstellar envelopes, 
from which we conclude that the SMA sources are very likely protostars in the Class evo- 
lutionary phase. We cannot advance more in our interpretation of the individual sources 
due to a lack of complementary high resolution observations at infrared and (sub)millimeter 
wavelengths of this region (namely at 850 //m). We would require such observations to build 
individual spectral energy distributions and model these to determine protostellar envelope 
and disk masses. We are however able to discuss the fragmentation of D-MM1 and the for- 
mation of these sources. There are 11 sources within D-MM1 (excluding one foreground star 



Young et al. 2006), 7 of which are the SMA sources we are reporting in this Letter. These 



sources are more densely packed than those in the Spokes cluster, which have a characteristic 
spacing similar to the Jeans length (Teixeira et al. 2006). Our interferometric observations 
recover ^90% of the flux of D-MM1 measured by Peretto et al. (2006) using their back- 
ground subtracted maps. This implies that the single dish observations correspond to the 



convolution of the individual fluxes emitted by the SMA sources (as suggested by Wolf-Chase 



et al. 2003). If we assume that each protostar has 0.5 M , then a rough estimate of the pre- 



fragmentation density of D-MM1 can be calculated by adding up the measured masses of 
the envelopes and the adopted protostellar masses, which gives us 10.9 M , and distributing 
this total mass within a region of 7870 AU radius (Peretto et al. 2006). The density we 
estimate in this manner is 8xl0 5 cm~ 3 . The Jeans length, Aj = a/ (irk^T) / (Gn) / (/imn) (A;b 
is the Boltzmann constant and G is the gravitation constant) we obtain for D-MM1 is 5.9" 
(0.023 pc or 4740 AU). Here we use a temperature T of 10 K since it is logical to assume that 
the pre-fragmented starless core would be colder than the present temperature of D-MM1, 
23 K. The corresponding Jeans length is very similar to the mean distance between the 11 
sources, 6.9" (0.026 pc or 5520 AU). However, the mean nearest neighbor separation between 
the sources is 2.3" (0.009 pc or 1840 AU), so it is not entirely clear if the fragmentation of 
D-MM1 was purely thermal. It is interesting to compare our results for D-MM1 with the 
northern subgroup of the Serpens cluster, which consists of a rich clustering of Class and 
Class I sources. Placing these Serpens sources at the distance of NGC2264 and scaling the 
1.4 mm fluxes measured by Hogerheijde et al. (1999) accordingly gives us fluxes comparable 
to what we are measuring for our SMA sources. Winston (2007) found that the average sep- 
aration of the Serpens protostars is about 5000 AU with some of the sources being as close as 
2000 AU. They also calculate the Jeans length for the region to be 0.024 pc, which is similar 
(neglecting projection effects) to what we find. If this is a characteristic scale then this could 
be an indication that thermal physics is the underlying engine that regulates star formation, 
in at least some highly clustered environments. Given that the nearest neighbor separations 
of the sources in D-MM1 are smaller than typical envelope sizes (5000 to 10000 AU), these 
sources may be competing between themselves for the accretion of the surrounding material 
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of the core. Therefore, we cannot rule out the possibility that some kind of competitive 
accretion (Bonnell et al. 2001 Peretto et al. 2006) may be occurring in D-MM1. Additional 
data is required, such as the measurement of the velocity dispersion of the sources, before 
we are able to draw any conclusions regarding this issue. Our finding holds within it an 
interesting implication, specifically, that there are two scales of fragmentation in the Spokes 



cluster: one that formed D-MM1 and other bright members of the cluster (Teixeira et al. 



2006) and a second associated with the fragmentation of the D-MM1 core itself. These 



two fragmentation length scales are correlated with the mean density in those regions, as 
would be expected if fragmentation was dominated by thermal physics. The Jeans mass for 
D-MM1, Mj = 17T 3 / 2 n" 1 / 2 M , is 0.6 M . Since D-MM1 would have had 10.9 M in its 
pre-fragmented state, we expect the core to have fragmented into about 17 sources if ther- 
mal pressure was the dominant support against gravitational collapse. This value is within a 
factor of two of the number of sources in D-MM1, 11. A more detailed analysis is obviously 
required in determining the underlying stellar masses. It is however interesting to note that 
if the Class 0/1 sources in D-MM1 have in fact M~O.5M , then their bolometric luminosi- 
ties could be similar to those of IRAS 05173-0555, RN043MM, CB230, IRAS 18148-0440, 



or RN015FIR, i.e., 7-10 L (Froebrich 2005), which would yield a combined luminosity of 
77-110 L . The bolometric luminosity of D-MM1 is 107.5 L (|Wolf-Chase et al.|[2003|. 



4.2. X-ray emission associated with SMA-1 



Flaccomio et al. (2006) used the Advanced CCD Imaging Spectrometer (ACIS) onboard 
the Chandra X-ray Observatory (CXO) to obtain a 97 ks long exposure of NGC2264. We 
searched their X-ray catalog and found that source # 237 (CXO ANC J064105.5+093407: 
(a, 5)(J2000) = (06 h 41 m 05.55 s , +09°34 / 07.68") is located approximately 0.6" away from SMA- 
1. We note that the SMA has subarcsecond pointing accuracy and the CXO has a positional 
accuracy of typically 0.6". To our present knowledge there are no other infrared or millimeter 
sources positioned closer to CXO ANC J064105. 5+093407 so it could potentially be associ- 
ated with SMA-1. The ACIS data has been recently re- analyzed and made publicly available 
through AN archive of CHandra Observations of Regions of Star formation (ANCHORS) 



( |http: / /cxc.harvard.edu /ANCHORS f\ |Wolk fc Spitzt^|2006fr . ANCHORS is a web based 



archive of all the point sources observed during Chandra observations of regions of star for- 
mation and their fitted spectra. Although the number of photons detected (9) was very low 
it is nontheless illustrative to mention the spectral fit obtained to help gain some insight into 
the nature of the source. The fitted spectra of CXO ANC J064105. 5+093407 corresponds to 
an absorption hydrogen column density of Nh=16.58±5.21x10 22 cm -2 , kT=0.90±0.13 keV 
(i.e., a plasma temperature of 10.4+1.5 MK), and a 0.3 solar metal abundance. We de- 
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rive an equivalent visual extinction Ay of 103±32 magnitudes using the gas-to-dust ratio 
A^tcm" 2 ) = 1.6 x 10 21 A v (mag) ( |Vuong et alTj|2003| ). Dust extinction maps of NGC2264 
using the near-infrared color excess method (e.g. Teixeira et al.||2005 and references therein) 
show that the region of the Spokes cluster is highly extincted, with Ay > 40 magnitudes, 
unpublished data. So the measured Ay of CXO ANC J064105. 5+093407 is consistent with 
it being either an Active Galaxy Nucleus (AGN) observed through the dense molecular cloud 
or a deeply embedded young stellar object. The fitted spectra has very little emission in the 
soft X-ray energy regime (0.5-2 keV) and is mostly hard X-ray (2-10 keV), peaking at 3keV. 
AGN spectra are characterized by hard X-ray emission (e.g. Mushotzky et al. 2000), how- 
ever, the plasma temperature derived from the fit is below that which is expected for AGNs 
(Tpi asma AGN ^ 10 9 K) and very similar to what is observed towards protostellar objects in the 
Coronet cluster [Forbrich et al.|fl200~6~] ). Therefore CXO ANC J064105.5+093407 is associated 
with a protostellar object, very likely SMA-1. Winston (2007) also reports Class 0/1 sources 
with X-ray emission. Unfortunately the CXO ACIS-I camera was pointed and orientated 
such that D-MM1 fell precisely in the gap between the chips of the array: this may be the 
reason why no other X-ray source is detected in D-MM1. 



5. Summary 

We discovered a dense micro-cluster of Class sources within D-MM1, with envelope 
masses ranging from 0.4 to 1.2 M and average radius of 600 AU. Five of the sources are 
associated with NIR/MIR emission, and one of the SMA sources, SMA-1, is associated with 
X-ray emission. The separations of the sources indicate that the fragmentation length scale 
of this core is significantly smaller than that of the Spokes cluster and comparable to the 
Jeans length in D-MM1. This is consistent with hierarchical thermal fragementation: with 
a primary fragmentation along the lower density (10 4 cm -3 ) filaments of the Spokes cluster, 
and a secondary fragmention of the higher density core D-MM1 (~10 6 cm -3 ). 



The SMA is a joint project between the SAO and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia 
Sinica. P. S. T. acknowledges the scholarship SFRH/BD/13984/2003 from the Fundagao para 
a Ciencia e Tecnologia (Portugal). 
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Fig. 1 — K s -band (a) and Spitzer IRAC 8 /an (b) images QYoung et al.|2006[ ) of D-MM1. The 
images are centered on (a, <5)(J2000) = (06 h 41 m 05.8 s , +09°34'08.8"). The SMA continuum 
data is overlayed in black contours ranging from 5 to 9 sigma (3.5 mJy beam -1 ) in both 
panels. The ellipse in the lower right corner shows the FWHP synthesized beam. 
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Table I. Properties of the SMA sources 



Measured Calculated 



ID. 


R.A» 


DEC>) 


SizeW 


P.A. 


F230GHz peak 


F230GHZ int. 


NIR 


R>gcom 


Mass 


n 




06 h 41 m 


09° 34' 


(") 


(°) 


(mJy beam -1 ) 


(mJy) 


/mir( c ) 


(AU) 


(Mq) 


(10 7 cm- 3 ) 


1 


05.60 s 


08.0" 


1.8(±0.6)x0.7 d 


70±13 


21±6 


29±14 


Nd, Mp 


448 


0.4±0.2 


14.7 


2 


05.64 s 


05.7" 


3.1(±0.7)x0.7 d 


23±6 


24±6 


56±19 


Nd, Md 


592 


0.7±0.2 


12.2 


3 


05.69 s 


06.7" 


2.6(±0.6)x0.9(±0.6) 


42±13 


27±6 


69±21 


Np, Mp 


608 


0.9±0.3 


13.8 


4 


05.87 s 


11.7" 


1.8(±1.0)xl.2(±1.2) 


96±45 


19±6 


41±19 


none 


584 


0.5±0.2 


9.2 


5 


05.98 s 


09.8" 


2.6(±0.6)xl.6(±0.5) 


59±32 


27±6 


90±25 


none 


816 


1.1±0.3 


7.5 


6 


06.09 s 


08.6" 


2.5(±0.5)xl.5(±0.4) 


87±18 


30±6 


94±24 


Md 


776 


1.2±0.3 


9.0 


7 


06.16 s 


09.5" 


2.1(±0.6)xl.2(±0.6) 


10±57 


21±6 


51±20 


Md 


632 


0.6±0.2 


9.1 



( a >J2000. 

W deconvolved FWHM size derived from fitting a 2D elliptical Gaussian to the continuum SMA sources. 

( c )Nd(Md): associated with diffuse K s -band(8 /im) emission; Np(Mp): associated with pointlike X s -band(8 /^m) emission. 

Wthe minor axis is unresolved and we use an upper limit of beamsize/2. 



